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Q Abstract 

The free energy due to the vacuum fluctuations of matter fields 
on a classical gravitational background is discussed, ft is shown 
explicitly how this energy is calculated for a non-minimally coupled 
scalar field in an arbitrary gravitational background, using the 
heat kernel method. The treatment of (self-)interacting fields of 
higher spin is outlined, using a meanfield approximation to the 
gaugefield when treating the gauge boson self interaction and the 
fermion-gauge boson interaction. 

1 Introduction 

When I first started, together with Prank Antonsen also of the NBI, on 
the study of the Casimir effect in curved space, I did so utterly ignorant 
about previous theoretical devellopments in the subject, thus our ap- 
proach came to be different from the standard one. The motivation for 
the study was our interest in wormholes and time machines and the fact 
that a wormhole almost has the shape of a cylinder and a cylinder looks 
as an obvious candidate for a calculation of the Casimir effect. Thus 
the original idea was to study Casimir driven space-time evolution (to 
see if wormholes are stable) and because the Casimir energy is finite, we 

^For a printed version of this paper, see the Proceedings of the 3rd International 
Alexander Friedman Seminar on Gravitation and Cosmology, Friedman Lab. Press, 
St. Petersburg 1995. 



1 



figured out that it had to be incoorporated in the Einstein equations (an 
insight, as prof. Grib kindly pointed out to us afterwards, other people 
had while wc were still in kinder-garten) . This paper will describe the 
method wc developed for determining the Casimir effect, taking as the 
starting point the flat space approach. 

With the action given by 

^ = -^Einstein-Hilbert + '^'matter (1) 
one obtains the classical Einstein equation 

-J^-G,.= T,. ^ (2) 

which, upon treating matter fields quantum mechanically gives follow- 
ing, 'corrected' (or first quantized) Einstein equations 



— ~ (JTILV =< J-UV >~ -T' (j) 



To get the right hand side of the first quantized Einstein equation in 
the case of a vacuous space-time one needs to determine the Casimir 
energy momentum tensor which in turn can be found from the Casimir 
free energy. Proceeding as in fiat space one expects to renormalize this 
quantity as follows 

-^ren — F -^Minkowski (4) 

where one subtracts the infinite free energy of the Minkowski vacuum. 
This prescribtion for regularization however, only works in locally flat 
space-times such as the hyperspatial tube. Prom a practical point of view 
there is no need to worry, though. By choosing the method described 
below one gets expressions that are readily zeta-renormalizable 

2 Determining the free energy from the zeta- 
function 

The Helmholz free energy is connected to the generating functional by 
the relationship 

F = -]-hiZ (5) 
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where P = (ksT)'^ is the inverse temperature and where the partition 
function is given by the functional integral which, in the case of a free 
scalar field, is 

Z = J e'^V^ = J e^/ v^'^'^5<^(°-'"')<^P0 = (det(D - m^))"^ (6) 

with the curved space d'Alembertian is given by 

□ ^ -1=5^(^5^-5.) (7) 

where g denotes the determinant of the metric. The use of a curved 
space d'Alembertian is the only change from the flat space approach. 

The determinant of an operator. A, is, by definition, the product of its 
eigenvalues 

det(^) = ]jA (8) 

A 

so with the zeta-function given by 

A 

these two quantities are related by the following equation 



ds 



=0 = ^ - InA • \-'\s=o = - ^In A = InJJ A = - lndet(y4) 
A A A 

(10) 



In the case of a minimally coupled scalar field, A = n — m^, one thus 
gets the following expression for the free energy 

3 Determining the zeta-function from the heat 
kernel 

Call the operator of interest A, which generally varies in space A = A{x), 
and the corresponding eigenfunctions ipx so that 

A^Px = AV'A (12) 
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The heat kernel, Ga{x,x' ,a), is the function satisfying the heat kernel 
equation 

d 

AGA{x,x';a) = -—GA{x,x';a) 
oa 

subject to the boundary condition Ga{x, x'; 0) = 6{x — x'). 
Note that this equation is satisfied by 



(13) 



Using this spectral representation of the heat kernel in turn makes it 
possible to show the following relationship between Ga and (a 



poo f rOO f 

J daa"'^ J y^d'^xGA{x,x,a) = J daa^~^ J y^d*x^\'ipxix)\'^e 



2-Xcr 



,s-lg-Ao- 



-A<T 



roo 

A ^ 

roo 

= yX-' d{Xa){Xay-\ 

= (EA-^)r(.) 



= ns)CA{s) 

where the fact that the V'a's are eigenfunctions of a Hermitean operator, 
and thus can be chosen to be an orthonormal base, has been used to 
perform the space-time integral. This gives us 

Ca{s) = J daa^~^ J GA{x,x;a)^/^d'^x 

and so, finally, for the minimally coupled scalar case 



gd'^x\i)x{x){'' 



(14) 



^□-m2 = 4|5=oCD-m2(s) 
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gdxdslg^o d(^^^Tr^GD-m^{x,x;a) 



r(s) 



roo 

= / ^n-m'^ix) 
Jo 



(15) 
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where J^{x) is the free energy density. Note that the integral over x is 
taken along the diagonal x = x' . 

Thus to calculate the zero-point energy of a scalar field one has to deter- 
mine the scalar field operator (the d'Alcmbcrtian) in the relevant space- 
time, solve the corresponding heat kernel equation and subsequently 
calculate the zeta-function from which the generating functional and 
thus all relevant quantities (including the zero-point energy) can be cal- 
culated. 

4 Determining the heat kernel for a non- Minimally 
Coupled scalar field 

Knowing how to relate the free energy (through the zeta-function) to 
the heat kernel let us proceed to actually determine the heat kernel for 
the scalar field operator, □ + with R being the curvature scalar (one 
can of course use any other function, including a mass term; —w?). 
To this end, first rewrite the d'Alembertian in terms of the local vierbeins 
(i.e. the coordinate base of a freely falling observer, = where 
Greek indices refer to general coordinates while Latin indices refer to 
the local inertial frame (the metric is g-^y = ?7a6e^e|^ where r]ab is the 
Minkowski metric, g is the determinant of the metric and e the vierbein 
determinant, e = ^/—g)'- 



where Dq = rj dad}, is the flat space d'Alembertian the heat kernel of 
which is known to be 
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= -9^(er?«''e^e|;a,)+Ci? 



= -e-a„(er?«''e(:e,^eX)+ei? 
= le-5^(er7«^e^e^eIJ)a„ + r?"^. 



(16) 



Gq{x,x' ;a) = {Ana) 




(17) 
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where A{x,x') = ds { = x — x' in Cartesian coordinates) where 
X, x' refer to freely faUing coordinates which one eventually may have to 
express in terms of 'general' coordinates. Q 

Now in the heat kernel equation 

□o + ^e™(9™(ee^^))9'^ + ^R^ Gn+^^j(x, x'; a) = -d„Gn+i^R{x, x'; a) 

(18) 

one can remove the first order term by substituting 

Gn+i^R = G(x,x';^7)e-^/^M (^'"(-i;))^-" (19) 

to get the following equation (which is not quite a heat kernel equation 
because the solution must respect the boundary condition G{x,x';a) = 
5{x — x') of the original equation, not G{x, x'\ a) = 6{x — x')) 



1 /I V 1„„ /I 



G{x,x';a) ~ —daG{x,x';u) 
(20) 



written compactly as 

(□o + /(x))G(x, x'- a) = -d^G{x, x'; a) (21) 

To determine G{x,x';a) from this equation, substitute 

G{x, x';a) = Go{x, x'; a)e^(^'^''") (22) 

to obtain 

□oT + dpTdPT + f + 2dp In(Go) • dPT = -d^T (23) 

where the last term on the left hand side vanishes on the diagonal, 
X = x' , because there, cf equation (17), the flat space heat kernel, Go = 

^ If one needs explicitly to use Dq in a calculation one must translate it so that it 
operates on the curved coordinates of which everything else depend. The translation 
reads Qo = ri'^^dadb = fj'^^eadfiieldu) = g^'^d^di, -\- rf'^e'i_(d^ie'^)du. However, in the 
following, it is perfectly possible to stay in freely falling coordinates and use the 
corresponding heat kernel of Do. This hinges on the fact that we need only evaluate 
the heat kernel on the diagonal, x = x' . 
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Go{x,x;a) is constant. 

Then, Taylor expand T in powers of a; 

oo 

T=Y1 TnU^ (24) 

rt=0 

The coefficents, t„, are easily determined: 

Prom the boundary condition we demand that (cf equations (19,22) 

^□+eil(^,x';0) = Go(x,x';0)e^(-'-''°)e-^/^^(^'"(-"))'^-" =(5(x-xO 

(25) 

Adding the fact that the flat space heat kernel already satisfies the 
boundary condition by itself, Gq{x,x';G) = S{x — x'), gives 

ro = \J le^{dm{ee^J)dx- (26) 

and collecting the cr^-terms one obtains 

n = -/ - Doro - dpTodPro (27) 

All higher order coefficents can be described by the recursion formula 
one obtains by inserting equation (24) into equation (23), remembering 
that the last term on the left hand side vanishes, and collecting the a"' 
terms 

Tn+l = TT ( °OTn + Y] dpTn-n'd^Tn' ) ;n > 1 (28) 

yielding for the next two coefficents 

T2 = f + higher order terms (29) 

and ^ 

T3 = —-dpfd^f + higher order terms (30) 

and so, from equations (19,22,24,26,27,29,30)), we finally have 
(47r(7j^ 

= J^^e'^" (l - (°oro + dpTodPTo)a + ^Do/a^ - ^dpfdPfa"^^ 
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giving the following renormalizcd expression for the Casimir free energy 
of a scalar field in an arbitrary gravitational background 



One should note that in the expansion (22,24) and subsequently the free 
energy one gets terms of higher and higher order in the curvature. Thus 
the term —fc in the heat kernel is the contribution from the classical 
gravitational background while the other terms that have been written 
out explicitly can be thought of as first order quantum corrections to 
that background. As we know general relativity to be renormalizable 
only to the one loop level it would thus probably be meaningless to 
continue the expansion (much) further. 

5 Determining the heat kernel of (Spin 1) Non- 
abehan Gauge Bosons 

In order to be able to carry out computations for higher spins as well, we 
want to relate the cases of vector bosons and spin-^ fcrmions to that of 
the scalar field case and as the treatment of gauge bosons carries some 
resemblance to that of the non-minimally coupled scalar field this case 
will be considered first. 

When considering a Yang-Mills field in a curved space-time then, in or- 
der to obtain the field strength tensor the naive guess is to replace the 
derivatives of the Minkowski space field tensor with covariant derivatives 
which is not correct as this leads to a non-gauge covariant expression 
(eventhough, accidentally, it gives the right answer in the case of abelian 
fields) . Instead proceed by considering the full theory of Dirac fermions 
interacting minimally with the gauge fields as well as with the gravi- 
tational field. In order to preserve local gauge and Lorentz covariance 




167r2/3 



J V^d\ ln(/) + - (DoTo + dpTodPTo)f ln(/) 

+ (DoTo + dpTodPro)/ - ^(Do/) ln(/) - ^(5p/5^'/)/-i 




(31) 
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construct a covariant derivative of the form 

Dm = e^^(5^ + |<(x)Xp, + igA^^ix)Ta) (32) 

where is the vierbein (local base vectors), ujj^'^(x) is the spin connec- 
tion being the gravitational analogue of the gauge field A'^{x) and Xpg 
the corresponding Lorentz group {SO (3,1)) generators analogeous to the 
gauge group generators T^. As before, greek indices refer to curvilinear 
coordinates while latin indices refer to local Lorentz coordinates and are 
also used for gauge group indices (it should be clear from context: in 
general small Latin letters from the begirming of the alphabet will be 
used to denote gauge indices, reserving letters from the last half of the 
alphabet for use as Lorentz indices). 

As in flat space field theory the gauge field tensor is obtained from 
the commutator of the covariant derivatives 

[Dm, D„] = Smn''{x)Dg + ^RmJt''{x)Xp, + iF^^Ta (33) 

(where S^ni^) is the torsion and i?^^*(x) is the Riemann curvature 
tensor). One obtains (after a lengthy calculation) the field strength 
tensor 

Kn = e^m<^{d,K - d.Al + igf^^A'^Al) (34) 

The generating functional of the full theory of a fermion interacting with 
a non-abelian gauge field is 

Z = J DA^ J Dip X)'0g'^BaM9e/ieW+'^/er-m.ion 

= j DA^j i;V'£'^e-i/^-^""^'^^"+'/^^'"^'"'^^'^"'(35) 

Refering back to equation (32) we see that the fermion part of the action 
contain reference to the gauge field so that one a priori cannot carry 
out the two integrations independently. One could probably treat the 
fermion action as a source term when doing the gauge field integration 
and then subsequently do the fermion integration. Instead however, 
I'm going to make a mean field approximation to A^ in the fermionic 
integral, as well as in the higher order terms of the bosonic integral 
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(see later) enabling one to consider the bosonic and the fermionic parts 
independently. 



The bosonic part of the generating functional, 

Z = J DA^e-y ^^-^^"^s""^' 



(36) 



can (using commutation relations, suitable normalization and the 
Lorentz condition) be given the form 

Z = JdA^ exp J ^d'xAl^ [-WdpdP + 5fe"(5„e'"'^ - d^e^)ePdp 

+gf,-^idnA"^^ - d^A'^^) + ic5'/e6c//e^;^^'] (37) 

In order to perform this path integral, make it Gaussian by choosing the 
following mean field approximation: 



DAa exp 



-gd^xAli- \-5t5^dpdP + 5t{dne^^^ - d^e^^)eld^ 



We now proceed to determine the mean fields. 

5.1 Determining the gauge mean field 

By definition^ 

I Al,{x)A\[x')e'Sj^A 



Al 



(38) 



{A^^{x)Ai{x')) 
It can be shown [1] that 

{A'i,{x)Ai{x')) 



J e'^VA 



(39) 



5^S 



-1 



2 \6A-^ix)5AUx') 



<5"Vn5(x-x')G(x,x') 
(40) 



^Where ever the mean value of a any odd power of the gauge fields occur, we 
substitute vC4^ for A. 
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where S denotes the appropriate action and where G{x, x') is a kind of 
Greens function related to the heat kernel G{x,x';a) by the relationship 

roo 

G{x,x') = - GA{x,x'-a)da (41) 
Jq 

The heat kernel is calculated using the same method as is done for the 
full theory in the following subsection so I just quote the lowest order 
result from [1] 

{A%,{x)Ai{x')) = (7 - m (-a^C^ + Isr^n^ + AT (42) 

where 7 is the Eulcr constant, S"^^ = — (9„e™^ — (?™e^)e^j and where 
is defined by equation (44) below. To explicitly find the meanfield start 
by putting it equal to zero in the action/formula (43) in order to calcu- 
late the first approximation to the meanfield. Then introduce this first 
approximation into the above formula to get the second approximation 
and so forth. 

5.2 The heat kernel for gauge bosons in the mean field 
approximation (cont.) 

The full heat kernel equation for the operator in the Gaussian path 
integral (39) is 



9_ 
4 



< IsTg^feycfr^lA^" >] Gt{^,x'-a) = -d,GZ{x,x' ■,a) 



2 

or, in short notation 

[5t,5'^dj>dP - 51, {dre^^ - a^e^e^Sp + /;r(< A >)] C^i^^ix, x'; a) = -d.Gl^{x, x 

(43) 

The first order term is eliminated, as before, by the substitution 

G = Ge^ /(9ne-''-a'"ej;)e^<ixp (44^ 

leading to the following equation 



9' 
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+ < \5^g''feVcfd\^lAP'' >] gL''(x,x» = -d,GZ{x.x';cT) 

which is of the form 

^ [5^5^, Do + aT] Gti^, x'- a) = -d,G^{x, x'; a) (45) 

The heat kernel becomes a matrix-valued function so assume G to be of 
the form 

(5?-(x,x» = G°(x,x'; la)(e"r(-'-''^))r (46) 

9 

where G°{x,x';a) denotes the heat-kernel of Dq (thus G°{x,x'; p-cr) is 

2 

the heat kernel of ^Qq) a-nd T is some matrix (T)^™ = Tf^^. Inserting 
this expression for the heat-kernel into the heat equation we arrive at 
an equation for T^^^"' 

□oT^r + {dpT§T) {&>tI') + or = - (47) 

where a summation over repeated indices is understood and where the 
fact that, along the diagonal x = x', the fiat space heat kernel is con- 
stant, and hence its derivatives vanish, has been used to eliminate the 
term 2dpG°dPT 

We will furthermore write as a Taylor series 

oo 

rr(^,^) = E^^tr(^K (48) 

i/=0 

Due to the boundary condition one has, from equations (45,47) and the 
fact that G° satisfy the boundary condition by itself, that 

d'ro r = -l^i {due"'' " ^'"O (49) 
and because of this and equation (48) the next coefHcent becomes 

\T = ^65p[(5ne™^-a-ej:)e^]-lj,-[(afce™'^-a"e0e^][(9„e'=^-a'^eit)e^^]+/r(< A 

(50) 



12 



Finally, inserting the expansion (49) into the equation (48) for T, yields 
a recursion relation for the higher order coefficients, t^'^^"'^, 

r^'\7 = -^,i^or^^^Z+ t id.r^''~'''^:rWr^''Vn) > 1 (51) 

6 Heat Kernel Equation and Zeta-Function for 
spin 1/2 fermions 

In order not to get calculational constipation, when relating the zeta- 
function of a fermionic field to that of a scalar one, I shall do the calcu- 
lation the Quick 'n' Dirty way (in equation (54)). 

6.1 Free spin 1/2 fermions 

The zeta-function of the operator A is 

a(.) = E^"' (52) 

A 

where A denotes the corresponding eigenvalues, Atpx = -^V'A) a-^d conse- 
quently, we also have the zeta-function for the operator A'^ (with eigen- 
values A^) 

Ca^{s)=^{XY' = Ca{2s) (53) 

A 

Now note that for a free fermion field the Dirac operator is 

y = 7™e^(5M + ^<(^)^P.) (54) 

Representing the SO{3, 1) generators in terms of the sigma matrices, 
'^PQ ~ i[7p'7g]' obtains for the derivative squared 

y = {D + ^fR) . U (55) 

(where I4 is the four dimensional unit matrix and = 1/8) establishing 
the link between the scalar and the fermion cases if one remembers 
to include a factor of 4 (one for each spinor component) in the Dirac 
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case (corresponding to taking the trace over the unit matrix, cf a smaU 
generaHsation of equation (8)). Thus we find 

Cy(^) = Cy^(^^)=4C&+«X|) (56) 

which, by use of section 4, makes it possible to determine the Casimir 
free energy of a free fermionic field if one remembers to take into account 
that for a Grassman field, the free energy is: 

F = -i(ln(det)0)^/2 (57) 

6.2 Fermions minimally coupled to a gauge field 

For real world purposes one should consider the case where the fermion 
field couples to a gauge boson field. The gauge invariant and space 
covariant derivative is 

= e^(5^ + + i5^^(x)T„) (58) 

one obtains, for the derivative squared; 

^2 = a + ^fR + 2gaP'iF^gTa + grfiA';,A''gTan + g{dPA;)Ta 

-Li^T ( rr i,,P1p^ A°- -I- 9/,)™"('p'^ A'^ — p^^ A'^ ) — 4m"'JpM A°- rr"^\ 

= U + ifR + 2gaP'iF;^Ta + gr,^'iAlA%n + g{A) (59) 
One can show that putting (the " carpet gauge" ) 

G{A) = (60) 

is an allowed gauge condition (by demonstrating that det(^^) ^ 0) [1]. 
Using the mean field approximation described in subsection 5.1, the 
gauge field dependent terms simply becomes a function which is in prin- 
ciple no diff'erent from the ^fR term to which it, for calculational pur- 
poses, can be added. 



14 



6.3 Fermions minimally coupled to a gauge field in the 
presence of (background) torsion 

The presence of (background) torsion, defined by equation (43), results 
in a small change in the calculation of the derivative squared of the 
preceeding subsection. Explicitly: 

^2 = ^^Dm^Dr, = V^'^DmDr, + ia^^ [Dm, + 7'"(^m7")i^n (61) 

which upon use of equation (43) and calculations as in the preceeding 
subsection becomes 

IJ)^ = U+ifR + 2gXXaP'^F;gTa + grf^ A^A'^Ta^ + G{A) + ia^^S^^^Dj, 

(62) 

Prom a calculational viewpoint the complication due to torsion thus 
consists of adding to the operator in the heat kernel equation a first 
order term which can be removed by the same token as it was done 
for the vector bosons (the presence of a sigma-matrix ensures the same 
complications) and two zero-order terms, one of them with a reference 
to the (meanfield of the) vector bosons. 

7 Conclusion 

I have outlined a method for explicitly calculating the Casimir energy 
for quantum fields on an arbitrary curved space background. The ex- 
pansion of the heat kernel and consequently of the zeta-function and 
Casimir energy involves higher and higher order derivatives of vicrbcins 
i.e. curvature (and of gauge fields, if present) which suggests that only 
the first few terms need be taken into account, in accordance with the 
results presented. In the case of interacting quantum fields I resorted to 
a meanfield approximation in order to be able perform the path integrals 
involved and a procedure for determining the meanfields was outlined. 
The first approximation to the mean field treats the quantum field as 
a free one, subsequent approximations do not. One would think this 
to be in line with what one would get calculating Feynman diagrams 
in curved space (could it be done) because interaction would involve a 
calculation to higher loop order and we therefore expect the results to 
be fairly reliable. 
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